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^ ! Abstract 
> 

Under natural assumptions a Feller type diffusion approximation is derived for criti- 
cal multi-type branching processes with immigration when the offspring mean matrix is 
primitive (in other words, positively regular). Namely, it is proved that a sequence of 
appropriately scaled random step functions formed from a sequence of critical primitive 



> 



OO I multi-type branching processes with immigration converges weakly towards a squared 

' Bessel process supported by a ray determined by the Perron vector of the offspring mean 

in 
o 

(N 



X 



matrix. 



1 Introduction 

Branching processes have a number of applications in biology, finance, economics, queueing 
theory etc., see e.g. Haccou, Jagers and Vatutin [5]. Many aspects of applications in epidemi- 
ology, genetics and cell kinetics were presented at the 2009 Badajoz Workshop on Branching 
Processes, see 



Let {Xk)k(zz^ be a single- type Galton- Watson branching process with immigration and 
with initial value Xq = 0. Suppose that it is critial, i.e., the offspring mean equals 1. 
Wei and Winnicki [21] proved a functional limit theorem A"^") X as n — )■ oo, where 
A'/"'' := n~^X\nt\ for t € R+, n G N, where \_x\ denotes the integer part of a; G M, and 
('^t)tGM+ is a (nonnegative) diffusion process with initial value Xq = Q and with generator 

(1.1) Lf\x) = mj'{x) + Iv^xfix), f G Cr(K+), 

where is the immigration mean, is the offspring variance, and C^(]R+) denotes the 
space of infinitely different iable functions on IR+ with compact support. The process (A:'t)feM+ 
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can also be characterized as the unique strong solution of the stochastic differential equation 
(SDE) 

dXt = m, dt + ^Jv^X+d'Wt, t E M+, 

with initial value Xq = 0, where (yVt)tm+ is a standard Wiener process, and x"*" denotes 
the positive part of x G M. Note that this so-called square-root process is also known as 
Cox-Ingersoll-Ross model in financial mathematics (see Musiela and Rutkowski [T71 p. 290]). 
In fact, {4:V^^ Xt)t£R+ is the square of a 4V^~^me-dimensional Bessel process started at (see 
Revuz and Yor [IHl XI.l.l]). 

Moreover, for critical Galton-Watson branching processes without immigration. Feller |1] 
proved the following diffusion approximation (see also Ethier and Kurtz [31 Theorem 9.1.3]). 
Consider a sequence of critical Galton-Watson branching processes (X^"^)^^^^, n G N, with- 
out immigration, with the same offspring distribution, and with initial value ^o"^ independent 
of the offspring variables such that u'^Xq^'^ — ^ ^ as ri — oo. Then X^""^ X as n ^ oo, 
where A'/"^ := n'^xf^Jj for t G M+, ri G N, and {Xt)tm+ 

is a (nonnegative) diffusion 
process with initial distribution fi and with generator given by (11. ip with = 0. 

A multi-type branching process {Xk)k&+ is referred to respectively as subcritical, critical 
or supercritical if g{m^) < 1, ^(m^) = 1 or g{m^) > 1, where Qim^) denotes the 
spectral radius of the offspring mean matrix (see, e.g., Athreya and Ney [1] or Quine 

[T8]). Joffe and Metivier [lOl Theorem 4.3.1] studied a sequence (X^"^)fcgz+ of critical multi- 
type branching processes with the same offspring distributions but without immigration if the 
offspring mean matrix is primitive and n'^X^^^ /x as n — )■ oo. They determined 
the limiting behavior of the martingale part (A4^"^)„gN given by Alj'^'' := ^fc=i -^IT'' 
with M'j^^ := X^^^ - E(x[^"^ | . . . , X^^}^) (see ([MD). Joffe and Metivier [TOl Theorem 

4.2.2] also studied a sequence {X^j^^)kez+, n eN, of multi-type branching processes without 
immigration which is nearly critical of special type, namely, when the offspring mean matrices 
Tn^"\ n G N, satisfy m^"'' = Ip + n~^C + o(?t,~-'^) as n — oo, and they proved that the 
sequence (n~-'^x|"|j)(g]R^ converges towards a diffusion process. 

The aim of the present paper is to obtain a joint generalization of the above mentioned results 
for critical multi-type branching processes with immigration. We succeeded to determine the 
asymptotic behavior of a sequence of critical multi-type branching processes with immigration 
and with the same offspring and immigration distributions if the offspring mean matrix is 
primitive and ^"^Xq"^ as n — oo, where /x is concentrated on the ray M+ ■ Um^-, 

where is the Perron eigenvector of the offspring mean matrix (see Theorem [3?T|) . It 

turned out that the limiting diffusion process is always one-dimensional in the sense that for all 
t G M+, the distribution of Xt is also concentrated on the ray R_|_-it^^. In fact, Xt = XtUm^, 
t G M+, where {Xt)t£R+ is again a squared Bessel process which is a continuous time and 
continuous state branching process with immigration. In the single-type case, Li [15] proved a 
result on the convergence of a sequence of discrete branching processes with immigration to a 
continuous branching process with immigration using appropriate time scaling which is different 
from our scaling. Later, Ma IT6l extended Li's result for two-tvue branching processes. Thev 



proved the convergence of the sequence of infinitesimal generators of single(two)-type branching 
processes with immigration towards the generator of the hmiting diffusion process which is a 
well-known technique in case of time- homogeneous Markov processes, see, e.g., Ethier and 
Kurtz |3]. Contrarily, our approach is based on the martingale method. It is interesting to note 
that Kesten and Stigum |T3] considered a supercritical multi-type branching process without 
immigration, with a fixed initial distribution and with primitive offspring mean matrix, and 
they proved that g{m^)~"'Xn — ?■ W almost surely as tt, — )■ oo, where the random vector W 
is also concentrated on the ray M+ ■ (see also Kurtz, Lyons, Pemantle and Peres [H]). 



2 Multi-type branching processes with immigration 



Let Z+, N, M, M+ and M++ denote the set of non- negative integers, positive integers, real 
numbers, non-negative real numbers and positive real numbers, respectively. Every random 
variable will be defined on a fixed probability space A, P). 

We will investigate a sequence {X^^^) f^^^^, n G N, of critical p-type branching processes 
with immigration sharing the same offspring and immigration distributions, but having possibly 
different initial distributions. For each n G N, /c G and i G {1, . . . ,p}, the number of 
individuals of type i in the A;*^ generation of the n^^ process is denoted by xj^^ . By 
^k^jii denote the number of type £ offspring produced by the j*^ individual who is of 
type i belonging to the {k — ly^ generation of the n}^ process. The number of type i 
immigrants in the k^^ generation of the n}^ process will be denoted by Consider the 

random vectors 



X 



X 



in)- 
kA 



An) 
>>k,j,i 



'Jn) 
^k,j,i,l 



An) 
_^k,j,i,p_ 



Sn) 



Sn)- 
'-kA 



-in) 
'k,p_ 



Then, for n, A; G N, we have 



(2.1) 



i=i j=i 



Here |Xq"\ '■ k,j G N, i G are supposed to be independent for all 

n G N. Moreover, {^fc"] j : k,j,n G N} for each 2 G {1, . . . ,p}, and {^^"^ : A;,?t, G N} are 
supposed to consist of identically distributed vectors. 

We suppose E(||^i"'^i jP) < oo for all i G {1, . . . ,p} and E(||£^"'^^ |p) < oo. Introduce the 



notations 



Vi. :=Var(cW,) G 



ppxp 



^pxp 



:= Var(£S^^) G 



me := E(£['^) G 



Note that many authors define the offspring mean matrix as mj. For k G Z+, let J-'^"'' 
a(xi"),xS"\...,xi")). By (EH), 



(2.2) E{xt^ I 4%) = Y: Xt\,m,^ + m, = m,xt\ + m.. 



i=l 

Consequently, 

(2.3) E(xi"^) = E(xt\) + m,, k,neN, 
which implies 

fc-i 

(2.4) E(X^"^) = m^E(x[,"^) +^m^me, k,neN. 

j=0 

Hence, the offspring mean matrix plays a crucial role in the asymptotic behavior of the 

sequence 

In what follows we recall some known facts about primitive nonnegative matrices. A matrix 
A e M^""^ is called primitive if there exists m e N such that A"" G M^+^. A matrix A e M!^""^ 
is primitive if and only if it is irreducible and has only one eigenvalue of maximum modulus; 
see, e.g., Horn and Johnson P, Definition 8.5.0, Theorem 8.5.2]. If a matrix A G MJ^^ is 
primitive then, by the Frobenius-Perron theorem (see, e.g., Horn and Johnson P, Theorems 
8.2.11 and 8.5.1]), the following assertions hold: 



• g{A) G g{A) is an eigenvalue of A, the algebraic and geometric multiplicities of 
g{A) equal 1 and the absolute values of the other eigenvalues of A are less than g{A). 

• Corresponding to the eigenvalue g{A) there exists a unique (right) eigenvector ua G 

called Perron vector, such that the sum of the coordinates of ua is 1. 

• Further, 

^(A)^"A" ^ := uav\ G M^+^ as n ^ oo, 

where va € ]R+_,_ is the unique left eigenvector corresponding to the eigenvalue q{A) 
with u\va = 1- 

• Moreover, there exist c^,r^ G with < 1 such that 
(2.5) ||^(A)-"A" - n^ll ^c^r^ for all n G N, 

where ||-B|| denotes the operator norm of a matrix B G R^^^ defined by ||-B|| : = 
sup||^ll=i \\Bx\\. 



A multi-type branching process with immigration will be called primitive if its offspring mean 
matrix me is primitive. Note that many authors call it positively regular. 



3 Convergence results 



A function / : R+ — MP is called cddldg if it is right continuous with left limits. Let 
D(M+,MP) and C(M+,]R^) denote the space of all M^-valued cadlag and continuous functions 
on M+, respectively. Let Poo(IR+,IR^) denote the Borel cr-algebra in D(M4.,]Rp) for the 
metric defined in Jacod and Shiryaev [9l Chapter VI, (1.26)] (with this metric D(]R_|.,Mp) is 
a complete and separable metric space). For M^- valued stochastic processes {yt)teM.+ and 
(3^i"^)teK+, n & N, with cadlag paths we write — ^ y if the distribution of 3^^"^ on 
the space (D(M+, M^), r'oo(K+, K^)) converges weakly to the distribution of 3^ on the space 
(D(M+,MP),r'oo(K+,MP)) as n -> oo. 

For each n G N, consider the random step processes 

Ar;"):=n-^x[2j, tGR+, n G N. 

For a vector a = (Q;j)i=i,...,p G M+, we will use notation a := J2^=i^i'^ti ^ MP^p, 
which is a positive semi-definite matrix, a mixture of the variance matrices V^j, . . . , V^p- 

Theorem 3.1 Let (X^""^)^^^ , G be a sequence of critical primitive p-type branching 
processes with immigration sharing the same offspring and immigration distributions, but having 
possibly different initial distributions, such that n~^xlj^^ — ^ XQUm^, where Xq is a nonneg- 
ative random variable with distribution /i. Suppose E(||Xq"''|P) = O(n^), E(||^^"'^| -p) < oo 
for all ie{l,...,p} and E(||£^^^ |p) < oo. Then 

(3.1) A-^") Xurr,^ as n ^ oo, 

where (A:'t)tei8+ is the unique weak solution (in the sense of probability law) of the SDE 

(3.2) dXt = vj^^m, dt + ^vl^^iu^^ V^)v^^X+dWt, t G M+, 
with initial distribution fx, where {'Wt)teM.+ is a standard Wiener process. 

Remark 1 We will carry out the proof of Theorem 13. II under the assumptions E(||^[^] .||^) < 
oo for all z G {1, . . . ,p} and E(||£[^-'||^) < oo. In fact, these higher moment assumptions are 
needed only for facilitating of checking the conditional Lindeberg condition, namely, condition 
(ii) of Theorem IA.3I for proving convergence (13. 4p of the martingale part. One can check the 
conditional Lindeberg condition under the weaker moment assumptions of Theorem 13.11 by the 
method of Ispany and Pap [7], see also this method in Barczy et al. [2]. If (i ^ 2 then it is not 
clear if one might get rid of the assumption E(|| Xq"^|P) = O(n^) in Theorem 13.11 

Remark 2 Under the assumptions of Theorem 13. 1| by the same method, one can also prove 
Xum^ as — )■ oo, where x\ := n'^i^X^^^.^ — m^"'*-'XQ"^), t G M+, n G N, and 
{^t)t&L+ is the unique strong solution of the SDE (13. 2p with initial value = 0. 



Remark 3 The SDE (13. 2p has a unique strong solution (A'/^°'*)tgiR^ for all initial values 
^■(1^°) = xoeR. Indeed, since \^/x — y/y\ ^ \/\x^-y\ for x,y ^0, the coefficient functions 



M 9 X I— )■ v^^m^ e M+ and M 3 x t-i- y f^^(w^ V^)vm^x~^ satisfy conditions of part (ii) 

of Theorem 3.5 in Chapter IX in Revuz and Yor [19] or the conditions of Proposition 5.2.13 
in Karatzas and Shreve [12]. Further, by the comparison theorem (see, e.g., Revuz and Yor 
[T9| Theorem 3.7, Chapter IX]), if the initial value A'q^°^ = xq is nonnegative, then A'/^^ is 
nonnegative for all t G M+ with probability one. Hence may be replaced by Xt under 
the square root in (13. 2p . 

Proof of Theorem \3 . 1[ In order to prove (13. ip . for each n G N, introduce the sequence 

(3.3) Mi:^) := Xi") - E(xi:^) | J^".)) = X^:^) - m,xtl, - m., k G N, 

which is a sequence of martingale differences with respect to the ffitration (-^i"^) ^g^^ • Consider 
the random step processes 

AlS"^=n-Mx;,") + ^M["M, tGM+, nGN. 
^ fc=i ^ 

First we will verify convergence 

(3.4) M^") Am as n ^ oo, 
where {M,t)teM+ is the unique weak solution of the SDE 



(3.5) dMi = Y^(n^,(Mt + tme))+0y^d>Vt, t G M+, 

with initial distribution jji := XQUm^, where {yVt)t!m+ is a standard p-dimensional Wiener 
process, aj"*" denotes the positive part of a; G and for a positive semi-definite matrix 
A G W^'P, \fA denotes its unique symmetric positive semi-definite square root. 



From (13. 3 P we obtain the recursion 

(3.6) Xf^ = m^x["\ + M^"^ + m„ ke N, 
implying 

k 

(3.7) X^^^ = mlxl^^ + J2 m\-\Mf^ + m,), k G N. 

i=i 

Applying a version of the continuous mapping theorem (see Appendix) together with (13. 4p and 
(13. 7p , in Section H] we show that 



(3.8) 



^(n) ^ as n ^ oo. 



where Xt := Ilm^{M.t+tm^) , t G IR+. Using 11^^ = Um^v'^^ and vJ^^Um^ = 1 we get that 
the process yt := v^^^Xt, t e R+, satisfies yt = vJ^^Ilm^{M.t + tm^) = vJ^^{M.t + ttrie), 
t G M+, hence Xt = ytUm^. By Ito's formula we obtain that {yt)teR+ satisfies the SDE (13.21) 
(see the analysis of the process (Vt'^"^)tm+ iii the first equation of (14. ip and in equation (14.21) ) 
such that 3^0 = '^m^'^o = '"m^i'^o'u-m^) = '^0: thus we conclude the statement of Theorem 
O ' □ 

Remark 4 By Ito's formula, the limit process {Xt)tm+ iii (13. ip can also be characterized as 
a weak solution of the SDE 



(3.9) dXt = nm,m,dt + nm,^JXt QV^dWt, teR+, 

with initial distribution II^^A^o = Um^{XoUm^) = Xou^^, since 11^^ n^^ = Um^v^u^ 

Remark 5 The generator of {AA.t)tm+ is given by 

Uf{x) = \{mrr.,{x + tm,)) V^]V, V)f{x) 



i=i j=i 

where V^^ij := (Cov((^i i i ^j))^=i ...^^ G M^. (Joffe and Metivier [TOl Theorem 4.3.1] also 
obtained this generator with = deriving (13.41) for processes without immigration.) 

4 Proof of M^""^ M and Af^^^ X as n ^ cxd 

First we prove 

Jl^ m applying Theorem for U = M, C/J,"^ = n-^xj,"^ and 
LtJ^") = n-^M'j^^ for n, /c G N, and with coefficient function 7 : M+ x ^ Rp^p of the 



SDE dSSD given by 7(^5 — y (nj^^(a; + tTn,^))~^ © 1^^. The aim of the following discussion 
is to show that the SDE (13. 5p has a unique strong solution (•^i^°OteR+ "^i^h initial value 
A4[,^°^ = ^0 for all VqEW. First suppose that the SDE (13.50 . which can also be written in 
the form 



dMt = J{vl, {Mt + tm,))-^{u^^ © V^) dW 



has a strong solution {-^i^^^) t,zM.+ with A^q^"'' = t/q. Then, by Ito's formula, the process 
^pivo)^ Q[yo^)^^^^, defined by 

p(^o) := «;^^(A4l^°) + tm,), Q?") := Mi"'^ - P^^^^u^, 
is a strong solution of the SDE 



(4.1) 



dVt = vl^^m.dt + ^/Vtvl^^^u^^QV^dWt, 
da, = -n^.m.dt+ x/vf (I^ - n^-l./w^. « Vc dW. 



with initial value {Vq^°\ Qo'°^) = {'^m^yo^ i^d — ^m()yo), where Ip denotes the p- 



dimensional unit matrix. The SDE (14. ip has a unique strong solution {vi^'\ Q'f''^) 



with an arbitrary initial value [Vq , Qq ) = (PO)Qo) ^ 
04. ip can be written in the form 



X 



since the first equation of 



(4.2) 

with b := vZ, nrii, G 



dVt = &dt+ ^/P+dW^ 



and 



where {'Wt)tm+ is a standard one- dimensional Wiener process. (Equation ( 14. 2p can be dis- 
cussed as equation ([S2D in Remark (^.) If (Pf^"\ QS^''^)tGR+ is the unique strong solution of 
the SDE dO) with the initial value {V^^"\ QJ^"^) = {vJ^^Vo, {Ip - Um^)yo), then, again by 
Ito's formula, 



is a strong solution of (13. 5p with A4o^"^ = j/g- Consequently, (13.5p admits a unique strong 
solution (A^l^"^)igjj+ -^o^"^ = for all ?/o e M^. 

Now we show that conditions (i) and (ii) of Theorem IA.3I hold. We have to check that, for 
each T > 0, 



(4.3) 



(4.4) 



sup 

*e[o,T] 



lj:E[Mi")(Mi"))^|^S]- / {n^:^UdsQV. 
k=i -^0 

- Y: E(||Mi")|pl I 4t) ^ for all ^ > 



k=l 



as 72 — )■ cxD, where the process (7?.j"'')fgiR , is defined by 



(4.5) 

By (E3D, 



ni"^ := (Air^ + tm^) , t e M+, n G N. 



(n) 



^ ^ k=l 



where we used that 



n^ mt = ( lim rrit ) mt = lim m^~^^ = 11, 



m^) -|- tm^ 



implies n^^(/p - m^) = 0. Thus = and 



Using (IA.4p . we obtain 



k=l k=l 

Hence, in order to show (14. 3p . it suffices to prove 

II I' r _ TT ^\^(")|l Jil^n o„r. 

,2 II lnt\ 



(4.6) sup ^ V ||(/^-n^ Ao, sup 

as n ^ CO. Using (13. 7p and Ilm^m^ = Um^, we obtain 



Hence by fl23|l . 

[nij-l [ntJ-1 fc 

k=0 k=0 k=l j=l 

^-^^l\\X^:^\\ + [nt\.\\m^\\+ J2 WMf^W]- 

V j=i ' 

Moreover, by (ETl) and flXSl) . 



||x(")|| + LntJ.||m,||+5^||Mf)|| , 
V j=i / 

where (7,^^ is defined by ( lA.Sp . Consequently, in order to prove (14. 6p . it suffices to show 
4. y ||M;")|| a 0, All^o'^ll ^0 as n ^ oo. 

^2 II J II ^2 11 U II 

In fact, assumption rT^X.^^^ — ^ \i implies the second convergence, while Lemma [A. 2 1 yields 
s^VnT\ EfilM^^in ^ 0. thus we obtain (lOD. 



Next we check condition fl4.4p . We have 



E(l|Mrfl,|Mi-||>n.} 



Moreover, n"^ Ei=i E(||M^"^f ) ^ as n ^ oo, since E{\\M)^>\\^) = 0{{k + nf) by 
Lemma IA.2I Hence we obtain (14.4p . 

Now we turn to prove (!3l8|) applying Lemma |A31 By (ISTTj) . A'^"^ = \I/„(A1^'"^), where the 
mapping ^„ : D(R+,Mp) -> D(M+,Mp) is given by 



[ntj 



vl/„(/)(t) :=mfJ/(0) + 5^m, 



[nti-j 



n 



for / e D{R+,W), t e M+, n e N. Further, X = ^(Al), where the mapping 
D(M+, RP) D(R+, RP) is given by 



vI/(/)(t):=n^.(/(t) + tm 



/gD( 



.,RP), 



t e 



Measurabihty of the mappings n eN, and \I/ can be checked as in Barczy et ah [2]. 

The aim of the following discussion is to show that the set C := {/ G C(R+,Rp) : 
n^,/(0) = /(O)} satisfies C G V^{R+,W), C C and P(A4 g C) = L 

Note that / G C implies /(O) G R ■ Um^. 

First note that C = C(R+,Rp) n TiQ^{{Ip - n^^)~i({0})) , where ttq : D(R+,Rp) ^ Rp 
denotes the projection defined by 7ro(/) := /(O) for / G D(R+,R?'). Using that C(R+,Rp) G 
Poo (see, e.g., Ethier and Kurtz [21 Problem 3.11.25]), the mapping W 3 (/p — n^^)x G 
RP is measurable and that ttq is measurable (see, e.g., Ethier and Kurtz [3l Proposition 3.7.1]), 
we obtain C g r'oo(M+, R^). 

Fix a function f E C and a sequence {fn)neN in D(Rp) with /„ /. By the definition 
of \E', we have \1'(/) G C(R'''). Further, we can write 

*n(/n)(t) = n^, (/„ + ^m,) + (mf J - n^J/(o) 



[nt\ 



hence we have 

||M/„(/„)(t)-vi/(/)(t)|| ^ iin^^i 



J-1 



77, 



1 

n 



m 



n 



fn 



1 „ 

— 



m 



[nt\ 



n 



/n(0) 



fn 



J -I 



n 



n 



For all T > and t E [0,T], 



fn 



[nt\ 



n 



fit) 



fn 



[nt\ 



n 



f 



[nt\ 



n 



f 



\nt\ 



n 



fit) 



^ujTif,n-')+ sup \\f^{t) - f{t)\\, 

t£[0,T] 



where uxif,-) is the modulus of continuity of / on [0,T], and we have UTif,n. ^) — 
since / is continuous (see, e.g., Jacod and Shiryaev [Q] VI. 1.6]). In a similar way. 



fn 



fn 



n 



^ujTif,n~') + 2 sup 

t&[0,T] 



By 



Further, 



[nt] 

El 



m 



[nt\-j 



n 



[nT\ 



\ — r. 



\nt\ 



{rrit — 



I \nt\ 

smce yrrv-^ — 
lim^^oo rn^''^^ = 



n^J/n(o)|| ^ [[(mL'^'J -n^J(/„(o) - /(0))|| + ||(mL"*J -n^J/(0) 

^c^, sup 

t£[0,T] 

nmJ/(0) = for all t e R+. Indeed, m^H^^ = m^lim^ 

and /(0)=n^^/(0) imply mL"*J/(0) = mL"*Jn^^/(0) = n^^/(0). 
Thus we conclude C C C'*,('f„)„gN- 

By the definition of a weak solution (see, e.g., Jacod and Shiryaev [HI Definition 2.24, 
Chapter III]), has almost sure continuous sample paths, so we have P(A4 G C) = 1. 

Consequently, by Lemma EH we obtain AT*^''^ = ^„(M^"^) ^(M) = A' as n oo. □ 
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Appendix 

In the proof of Theorem 13.11 we will use some facts about the first and second order moments 
of the senuences fxl^-'lfccy. , and ( m{!'^)uck. 



Lemma A.l Under the assumptions of Theorem \3 . 1\ we have for all k,n & N 



k-l 



(A.l) E(xi")) = ml E(X[,")) + J2 ^i^e 

3=0 



k-l 



(A.2) 



j=0 



fc-2 k-j-2 

+ m|(Var(X;,")))(mJ)^ + ^m^^ ^ [(m|m,) V^] (mj^. 

j=0 f=0 



(A.3) E(Mi^") I =0 for k,ne N, 



(A4) E[Mr)(Mfy|^l:4,^_,] 



if k^e. 



Further, 

(A.5) E(Mi''^) = for keN, 

ye + E(xt\)0V^ ^/ fc = 

zf k^ 



(A.6) E[m[")(mJ"V] 



Proo/. We have already proved (^K^, see dH. The equahty m["^ = X^"^ - E(xt"^ | J-^) 
clearly implies fjA3l) and fjA5|) . By ([211]) and (JSJl), 



Y(n) 



(A.7) Mi") = Xl") - 5^ A^,, E(^«,J - = (s, - E(e,)) + E E (^S. " ^(^11))- 

i=l j=l j=l 

For each fc, n G N, the random vectors {^1"-^ -E(|2 J, e^^^ -E(£[")) : j G N, z G {1, . . .,p}} 

(n) 

are independent of each others, independent of ^^-ii ^^"^ have zero mean, thus in case 
k = t we conclude ^KM and hence (jAj]). If A; < £ then Y.[Mf\Mf'^y = 
M^") E[(Mf V I J^il\] = by (lA3!l . thus we obtain (lAill and (lA6l) in case k ^ i. 

By (13 .yp and flA.ip . we conclude 

Xin) _ E(X1")) = m|(x(") - E(X("))) + E^-^'^f- 



Now by flXej) . 

j=i i=i 

k 

Finally, using the expression in f lXej) for Yj[Mf\Mf'^y] we obtain (jD). □ 
Lemma A. 2 Under the assumptions of Theorern VJ . 1\ we have 

E(||xi") II) = 0{k + n), E(||xi-)f ) = 0{{k + n)% 
mMt\\) = 0{{k + nfl% Hmt^r) = OHk + ny. 

Proof. By flXTj) . 

II E(xi"))|| ^ ||m||| ■ EdlX^'^^ll) + J2 IIKll • ll^^ll ^ Cm,{VCn + \\m^\\k), 

j=0 

where 

(A.8) Cm^ := sup ||mp < oo, C := sup n"^ E(||X{,"^f ) < oo, 

since (12 .Sp implies Cm^^ ^ c^^j^ + ||njn^||. Hence, we obtain E(||x[,"'' ||) ^ p|| E(X^"'')|| = 
0{k + n). 

We have 

E(||Mt")||) ^ Ve(||m[")P) = ^E[tr(Mi")(Mi'^y)] = y^tr [V, + E(Xa) V^] 

^ V^MV:) + y^tr[E(xt\) Vg], 

hence we obtain E(||M^"^ ||) = 0((A; + n)i/2) from E(||X^"^||) = 0(A; + n). 
We have 

E(||xi")|p) = E[tr(xl")(Xr)n] =tr(Var(Xr)))+tr[E(Xr))E(^^^^ 



where tT[E{X'^^)E{X^^Y] = ||E(Xr^)f ^ [E{\\X^^^\\)Y = Oi{k + n)'). Moreover, 



tr(Var(X^"0) = 0{{k + nY). Indeed, by (jAj]) and flXSjl 



II Var(X 



<E (11^.11 + 11^ 



,l|-||m,^-^-^||.E(||X^' 



ill 2 



k-j-2 



+ ||Var(Xo")||-||m|f + ||m,||-||V^||^||m^^f ^ ||m| 



i=o 



i=0 



< (l|V,||+C^J|Vdl-E(||x(")||))c^^A: 



+ (E(||x(")f ) + [E(||Xo"1|)]^)C^ Jlm.ll ■ IIV^IIP, 



where := X]i'=i II^^JI' hence we obtain E(||X 



(")||2^ 



0{{k + nY). 



By (ED, 



Mf\\^\\sf-E{et^)\\+Y: 



1=1 



X 



(n) 



hence 

E(||M 
Here 



^{p+ If E{\\e? - Eie?)f) + (p + 1)^ E | 



E 



X 



in) 



E 



i=l 



k-l.i 



X 



(n) 



2\ 2- 



,£=1 \ j=l 



/X 



(n) 



where 



E 



A-n) 







)• 


T-(7l) 



with (E[(ei;,^, - Ei^iljf])' < mi'Le - E(eil,,))1, hence 



E 



/X 



(") 



^E[(e;i,-E(e};i^,))iE[(xn 



(1) ^^4l 



\2l 



Consequently, E(||x("^ll2 



0{{k + nf) implies EdlM^^^n^^ 



k II ;-v.vv-^'v ; -VII — II j-0{ik + nf). □ 

Next we recall a result about convergence of random step processes towards a diffusion 
process, see Ispany and Pap |H1 Corollary 2.2]. 

Theorem A. 3 Let 7 : M_|_ x — > MP^^ be a continuous function. Assume that uniqueness 
in the sense of probability law holds for the SDE 



(A.9) 



dUt = i{tMt)^y\^u te 



with initial value Uq = Uq for all uq G W\ where {Wt)te_R+ is an r- dimensional standard 
Wiener process. Let fx be a probability measure on (W',B{{M'^)), and let (Ut)tm+ be a 
solution of flA.9p with initial distribution fi. 

For each n G N, let {U^l^^)kez+ be a sequence of p- dimensional martingale differences 
with respect to a filtration {J^l^'^)k&+- Let 



lnt\ 



k=0 

Suppose E(||C/["^f) < 00 for all n,k & N, and C/g"'' /^. Suppose that, for each T > 0, 



[1) sup 

te[o,T] 



lnt\ 

EE 

k=l 



0, 



[nTJ 

(ii) E E(||C/i")fl, 

k=l 



^k-i) ^0 for all e>Q, 



where 



denotes convergence in probability. Then U 



(n) ^ 



U as n — !■ 00. 



Now we recall a version of the continuous mapping theorem. 

For functions / and /„, n G N, in D(M+,Mp), we write /n — > / if {fn)nm converges 
to / locally uniformly, i.e., if supjg[o,T] "~ /(^)ll n -> 00 for all T > 0. For 

measurable mappings $ : D(M+,Mp) ^ D(M+,M5) and : D(M+,Mp) D(M+,M'?), n G N, 
we will denote by C'$,(#„),^g^ the set of all functions / G C(M+, W) for which $n(/n) — ^ ^(/) 
whenever /„ ^ / with /„ G D(M+,Mp), n G N. 

Lemma A. 4 Let {Ut)tm+ o.nd {lA.f^)t<m+, n E'N, be W -valued stochastic processes with 
cadlag paths such that W^") U. Let $ : D(M+,Mp) D(M+,M'') and : D(M+,Mp) ^ 
D(M+,M''), n G N, be measurable mappings such that there exists C C C^^, with 
CgPoo(K+,Mp) and P(i^GC) = l. T/ien A 



Lemma [A. 41 can be considered as a consequence of Theorem 3.27 in Kallenberg [TT], and we 
note that a proof of this lemma can also be found in Isnanv and Pan l8l Lemma 3.11. 
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